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Abstract 

We show that contrary to the claim made by Hallin and Liljenberg in Phys. Rev. 
D52 1150,(1995), flhep-th/9412188| ) the thermal correction to the thermal decay 
■ or pair production rate for a system placed in a heat bath in the presence of an 

ON , external electric field, is always nonzero in the finite as well as infinite time limit. 

Using the formalism outlined there, we reestimate the decay rate for different values 
of temperature, mass and time. We also try to identify the parameter ranges where 
. the quantity of interest agrees with that computed previously, at high temperature 

(in the infinite time limit), from the imaginary part of the effective action using 
£NJ ■ imaginary time and real time formalism of thermal field theory. We also point out 

that in the strictly infinite time limit, the correct decay rate as obtained from the 
work of Hallin et. al. tends to diverge. 
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1 Introduction. 

o 

The study of pair-creation phenomena in the presence of an external field in a heat bath 
has attracted some attention in the last few years. There have been attempts to find 
out the temperature dependence of the decay rate by various groups (see for instance 
references [JXJ] to [[|). 

There exists many formalisms to calculate the pair-production rate for such a system. 
One of them is to calculate the effective action, the imaginary part of which is related to 
the vacuum persistence probability and the other is to compute the density matrix using 
functional Schrodinger representation and from thence calculate the production rate. In 
the past, there had been investigations where the imaginary part of the effective action 
was calculated using the real time formalism or Thermo Field Dynamics (TFD)0 and 
also the imaginary time formalism @] of finite temperature field theory. On the other 
hand J. Hallin and Per Liljenberg had derived the vacuum persistence probability for such 
a system employing the elegant approach of functional Schrodinger representation-in ref- 
erence . Unfortunately the results of these investigations do not see to agree with each 
othre. Some of the studies have found a finite thermal correction to the production ampli- 
tude whereas other studies have claimed to have not found any. In references @ and [[§] 
a thermal correction to the decay rate have been found and the results of the two almost 
match except for a probable typographical error. 

On the other hand the results of references , || and || are in contradiction to [[| and 
H|]. In this note we shall concentrate on the findings of reference ||. The authors of 
reference || claim that at a fixed temperature, in the infinite time limit (as is assumed 
in the effective action formalism), the thermal correction to pair production rate vanishes 
thereby contradicting the earlier results found using the finite temperature effective action 
formalism. 
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In this note, starting from eqn. [208] of reference - obtained by the functional 
Schrodinger representation (for the fermionic case), we will show that there is a finite 
correction to the pair production rate, even in the very large time limit at any nonzero 
temperature. This is in direct contradiction to the results obtained in ||. It must be 
noted here that though under some approximations this result reduces to the results ob- 
tained in reference @ and || and we shall discuss the implications of that shortly. 

To be little more precise, we find that in the limit of very small time and high temperature 
the results of reference matches with those of B and found in the high temperature 
limit. On the other hand in the studies of H ana |3| the production rate is calculated in 
the infinite time limit as performed in the strict thermodynamic sense. In this paper we 
will also point out that the production rate when properly computed, turns out to be a 
function of time and it tends to diverge once the infinite time limit is assumed. 



2 Calculation 

In this section we start from the expression of probability P(0, 0; tf), of finding the system 
with no fermion anti-fermion pairs in a state after a time t = tf, as is given in eqn. [208] 
of Hallin and Liljenberg ||, 



P(0,0;t/) = exp 
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For the benefit of the readers following reference || we re-introduce the notations as given 
there in [Bfl 



th = tanh 



u 9 = ^A + (p3 + t f y 



A = (p 1 ) 2 + (p 2 ) 2 + (p3) 2 + (rn)\ v = V(eE) 
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t f = tf(eE) 
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With these definitions eqn. Q can easily be cast in the following form 



(2) 



P(0,0;t f ) = exp 
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Since the quantity inside the logarithm in eqn.|| is less than one we can expand the 
logarithm to arrive at 



P(0,0;t f ) = exp 



C(P)- 
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The expression for vacuum persistence probability, in the Schwinger sense is obtained 
from the eqn. |J as 



W = lim 



t f — >oo 



lnP(0,0,tf) 

vt f ' 



(5) 



In || the infinite time limit was taken inside the integrand before carrying out the inte- 
gration to find out the explicit tf dependence of the resulting expression. Though in some 
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cases this procedure might give identical result with the case where the limit is taken 
afterwards, here it does not seem to be so. Therefore, in order to compute the production 
rate W, unlike the authors of [H], we shall evaluate the integral retaining the explicit in- 
dependence before dividing it by tf and finally taking the limit tf — > oo. In order to 
simplify the situation a bit, from now on, we will work for a while with the integral inside 
eqn.0]. It is given by, 



oo r 

2/3o; 9 



2 dAe~ wA / _ 2 dp 3 1 - e 2/w = / _ dAe'^ \ _ n df - \ dAe^ K \ n dp 3 



J-f.. 2 7^2 ,/_r.. 2 

(6) 
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The first integral here gives the zero temperature Schwinger pair production probability 
and since the controversy in the literature is about the finite temperature piece we will 
concentrate only on the second term in eqn.H. In the second integral the quantity Uq 9 



is given by Uq = yA + (p 3 + tf) 2 and it also appears in the exponent thus making the 

evaluation with exact tf dependence more difficult. Therefore to circumvent this difficulty 
we propose to rewrite the exponential using an integral transform of the following form 

□ a 

e ~*= « r e ~£*f. ( 7) 

Now if we identify 2/3 with a and the quantity uJ Q 9 as y/s, we can write 

-/W = J_ [°° -u[A + (pS+t f r]-!*du 



20T Jo u % 

With the aid of eqn.||, we obtain from the second integral (eqn. |J) upon performing 
the A integration 

i = 4= r d f r d 4 n dAe-^e-^?y£ 

y/7T J-tf JO U2 J™? 

= -^=e m M dp 3 — r — - . 9 

V7T J-tf JO «2 (7T + U) 

Hereafter one can interchange the order of the integrations and carry out the p 3 integration 
first. Upon performing this integration the resulting expression can be cast in the form of 
an error function with y/u x tf as its argument. Rewriting the error function as a series 
we obtain 

8 2 00 (iA 2k ~ 1 (-l) k+1 r°° dwft k -W 2 a & 2 

i - p ,_ e -™ 2 -« / — — e~ u 4u (io) 
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In order to write the pair production rate in a more compact form, we make the following 
change of variables u — > -z(f ) 2 in eqn.JTU] and use the definition of pair production rate 
as given in eqn.|| without taking the limit tf — > oo before hand. And the resulting 
expression is, 



/_ -\ 2k-2 , h i 1 
eE 2 2 °° ( *f/5) (-1) + roo dzz (2k-5)/2 

" /3 2 (2tt) 2 v^F ^ (2fc-l)(fc-l)! io § + £ y 



2 This follows from the definition of the Bessel function. 
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It can be seen from eqn.fll l] that for any nonzero value of tf the pair production rate, 
contrary to the claim made in [Q, is always nonzero. At zero temperature, i.e, when 
oo, eqn.pl] goes to zero because of the (/3) 2 sitting in the exponent 



ft 



And it conforms 

to our general expectation that as temperature tends to zero the temperature dependent 
correction to the rate should also vanish. On the other hand at a nonzero temperature, 
in the thermodynamic sense, as tf — > oo the rate becomes divergent. This is an unique 
feature that comes out of the treatment of reference |J. This is the main result of this 
paper. In the remaining portion of the paper we shall try to evaluate this expression nu- 
merically for different values of the parameters and write down the approximate analytic 
forms for each of these cases. 



3 Discussion and Results 

In this section we evaluate W for three different values of the parameters namely (a) 

< 1, m < 1 and /3 ~ 0(1) or little more, (b) if ^ 1 but /3tf <C 1 and m <C 1 and 

lastly (c) both f3 and fa are <C 1 but if ^> 1. In the last case we would vary if to show 
how W changes with its variation. 



We shall start with ) and then move over to the other two cases. In order to 

derive an approximate analytical form for W for this case we start from eqn.|J instead of 
eqn.|IT| as given in this text. Since in this case both tf and fa are assumed to be 1 and 



the parameter (3 is of the order of one or slightly more, the quantity fty[ra 2 + (p 3 + tf) 2 } 

is very small. Since the integral there varies from to oo we would like to approximate 
the quantity (u + n) in the denominator by u to arrive at 
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Now we can perform the u integration by using the formula for modified Bessel-function, 

. r i - r i 

Jo 00 u^e'^-t du = 2 £ 2 K v 2^ . And using this form of the modified Bessel func- 



tion eqn.JT^] now comes out to be 
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(13) 



For the range of parameters that we are interested in the argument of the modified Bessel- 
function tends to be very small. In this limit one can replace the modified Bessel-function 



by its approximate form, Ks[z) 

ft -™, 2 
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. Using this relation eqn.|| becomes 



dp 3 2 
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For m < 1 the exponential can be approximated by 
in this limit, (using eqn.f5j) turns out to be, 



W 



eET 2 



2 



2(3 J [fa 2 + (f + tff 



(14) 



1. Hence the pair production rate 



(15) 



We have also evaluated, numerically, the pair-production rate from eqn.[Il| for tf = 
10~ 3 , with fa = 10~ 2 and have plotted it in figure fl[]. It shows the behavior of W as 
one moves from low temperature to high temperature. It can be seen from the curve that 
the general trend qualitatively agrees with that given by eqn.[15|. As has already been 
mentioned in the beginning, apart from the numerical factors, this result agrees with that 
of reference [@] and [[(§]] obtained in the high temperature and infinite time limit from 
the imaginary part of the thermal effective action. 
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Figure 1: Production rate vs temperature for tf = 1.0 x 10 3 and m = 1.0 x 10 2 . 

Having discussed the case for £j <C 1 and m < 1, we will try to evaluate the pair 
production rate for if ^> 1 but (3if <C 1 and fn -C 1. In order to do that, it will be 
convenient if we proceed from eqn.Jll]]. The first important thing to notice here is, as 
(3tf <C 1, the dominant contribution to the quantity of our interest comes from the first 
term in the summation. Hence, for the dominant contribution, it is sufficient to retain 
the k = 1 term where W can be approximated by 

W « 2{eE) ( r Z -^e~^-i\ (16) 

Since the dominant contribution to the exponential comes from z ~ (^;), which for the 

parameter range of our interest turns out to be quite large, we will proceed as before 
by approximating the denominator by z and then carrying out the z integration using 
the modified Bessel function of fractional order. On using the exact form of the Bessel 
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function eqn.[16] approximately comes out to be 

(eE) 



W 



(2tt) 2 /?< 



(17) 



We have estimated W numerically for m = 10 6 and tf = 10 3 for various values of 
temperature ranging from 10 6 to 10 7 and have plotted them in figure 0. It can be seen 
that the shape of the curve is basically due to the function e - *-" 1 ^ as is given in eqn.JT? 
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Figure 2: Production rate vs temperature for tf = 1.0 x 10 3 , m = l.d x 10~ 6 

Lastly in figure |§ we have plotted W for (3 = 10~ 4 and fh = 10 -6 as a function of tf, 
ranging from 100 to 101 in arbitrary units. One can see from this curve that as time tf 
increases W also keeps on increasing. So in the strict thermodynamic limit as tf — > oo 
the W would become divergent. 



4 Conclusion 

Starting from eqn. [208] of ref.plj we have shown that for a fermionic system, contrary to 
the claim made there (i.e. in |yj, at any non-zero temperature there is a finite correction 
to Schwinger's pair production rate for any value of tf.lt is shown that the final expression 



(i.e. eqn.JTT|]) comes as a power series in f3tf. Though for small values of (3 and tf, this 
expression seems to be finite, but in the proper thermodynamic sense it is not. Since for 
realistic situations (3 can never be zero, on the otherhand tf in principle can approach 
infinity, hence in the infinite time limit, the quantity W tends to diverge. The reason 
behind this apparently contradictory results seems to depend on how one takes the limit 
tf — > oo.In our view one should first carry out all the integrations and write down the final 
result as a function of tf, before taking the limit tf — > oo. Though we have carried out 
this detailed exercise for fermions only but in principle this can be extended for Bosons too. 
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We would like to mention here that the same quantity obtained using the finite temper- 
ature effective action formalism ( where infinite time limit is implicitly assumed) using 
both Thermo Field Dynamics and imaginary time formalisms ( see e.g. and 0]) yields 
a nonzero finite result. In view of the previous results we find this corrected result of 
reference Hto be quite amazing and in our view this deserves further study. Particularly 
one should try to understand why the different formalisms are giving different answers 
and which one of them is correct. But in case the corrected result of reference || which is 
obtained through the elegant formalism of Functional Schrodinger Representation turns 
out to be right (though we have some reservatio ns about it), it can certainly be used 
to estimate the thermal contribution to the soft processes dominated production mecha- 
nisms in the pre-equilibrium Quark Gluon Plasma production phase of relativistic heavy 
ion collisions, at the SPS and RHIC. 
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Figure 3: Production rate vs £/ for Temperature= 1.0 x 10~ 4 , m = l.d x 10~ 6 
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